Yukawa Interactions and Supersymmetric Electroweak Baryogenesis 
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We analyze the quantum transport equations for supersymmetric electroweak baryogenesis in- 
cluding previously neglected bottom and tau Yukawa interactions and show that they imply the 
presence of a previously unrecognized dependence of the cosmic baryon asymmetry on the spectrum 
of third generation quark and lepton superpartners. For fixed values of the CP-violating phases 
in the supersymmetric theory, the baryon asymmetry can vary in both magnitude and sign as a 
result of the squark and slepton mass dependence. For light, right-handed top and bottom quark 
superpartners, the baryon number creation can be driven primarily by interactions involving third 
generation leptons and their superpartners. 



If the universe was matter-antimatter symmetric at the 
end of the inflationary epoch then the microphysics of 
the subsequently evolving cosmos must have dynamically 
generated the cosmologically observed baryon asymme- 
try (the baryon to entropy density ratio tib/s), 

8.36 x 1CT 11 < n B /s < 9.32 x 1GT U (95% C.L.) (1) 

0. The Standard Model (SM) of particle physics and 
standard big bang cosmological model contain all the nec- 
essary ingredients (Sakharov criteria [l[) for successful 
baryogenesis: baryon number violation, charge conjuga- 
tion (C) and charge conjugation-parity (CP) violation, 
and departures from thermal equilibrium. However, fun- 
damental symmetry tests and collider constraints indi- 
cate that generating the observed tib/s requires physics 
beyond the SM (see e.g. 0|). 

Among most widely studied viable possibilities is elec- 
troweak baryogenesis (EWB) which is testable with 
low-energy searches for permanent electric dipole mo- 
ments (EDMs) and high-energy studies at the Large 
Hadron Collider (LHC) [4|. In this scenario, electroweak 
symmetry-breaking (EWSB) - a cosmological transition 
in which SU (2) l is broken at temperature T ~ 100 GeV 
- proceeds via a strong, first order phase transition dur- 
ing which bubbles of broken electroweak symmetry nu- 
cleate and expand in a background of unbroken symme- 
try. Particle-antiparticle asymmetries generated by CP- 
violating interactions at the bubble wall induce a non- 
zero density of left-handed fcrmions, ni e ft, that diffuses 
into the unbroken background where baryon number vio- 
lating SU (2) l sphaleron (electroweak sphaleron) transi- 
tions convert it into baryon number. The expanding bub- 
bles capture the non- vanishing baryon number and freeze 
it in by quenching the sphalerons, leading to 7Jr 
in the bubble interior. The baryon number is propor- 
tional to nieft, which in turn depends on CP-violating in- 
teractions and chemical potential equilibrating reactions 
such as Yukawa and nonperturbative SU(3) C transitions 
(strong sphalerons). 

In this Letter, we reanalyze ni e ft in the context of the 
minimal supersymmetric Standard Model (MSSM), one 
of the best motivated conjectures of physics beyond the 



SM, and observe new features which have been missed in 
previous works (see e.g. 

HIEIZLUS and Refs. therein): 
(1) Yukawa interactions between bottom quarks, Higgs 
bosons, and their superpartners cannot be neglected in 
EWB, even if the ratio of the vacuum expectation values 
(vevs) of the two MSSM Higgs doublets, tan/3 = v u /vd, 
is mildly larger than unity (a parameter region favored 
by current experimental constraints). This typically re- 
sults in a qualitative change from the standard picture: 
the first two generations of quarks and squarks decouple 
in EWB if the first two generations receive CP asymme- 
try mostly through strong sphalerons, as is true in all 
well known scenarios. (2) the MSSM prediction for tib/s 
can vary in magnitude and sign as the masses of the third 
generation sfermions are varied, even when the dominant 
source of CP violation is proportional to a single phase 
with a fixed sign. (3) there exist parameter regions in 
which left-handed (LH) leptons drive the baryon num- 
ber production, unlike the traditional situations in which 
EWB proceeds mainly through the interactions of the LH 
quarks with electroweak sphalerons. This occurs in the 
large tan/3 region, in which r Yukawa interactions are 
significant. Unlike standard thermal leptogencsis scenar- 
ios (see e.g. 0]), this new scenario does not require the 
participation of a right handed neutrino sector. 

In what follows, we first present the computational 
framework and analytic intuition for our main results. 
We subsequently give the full numerical results and dis- 
cuss their implications. 

Framework and Analytic Intuition. The current density 
jp for each particle species p satisfies a quantum Boltz- 
mann equation (QBE) of the form d\jp = Sp P + S^, 
where Sp P and S^ are, respectively, CP-conserving and 
CP-violating source terms which depend on the MSSM 
interactions and chemical potentials. Sp P includes the 
terms that push the system toward chemical equilibrium, 
while S^ contains the effects of CP-violating interac- 
tions involving the phase transition bubble. We have 
developed numerical solutions to the full set of QBEs for 
all MSSM particle species chemical potentials, including 
contributions to Sp p from previously neglected Yukawa 
and triscalar interactions as well as "supergauge" inter- 
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actions involving gauginos, particles, and sparticles. 

The full numerical results, which we report at the end 
of this Letter, can be understood by considering an ana- 
lytic solution which typically is valid in the limit of large 
tan /3 and superequilibrium which we explain shortly. In 
this regime, there exists a hierarchy of time scales in 
the phase transition dynamics that implies a set of sim- 
ple relations between particle chemical potentials and 
densities: Tdiff, associated with the diffusion of particle 
densities ahead of the advancing bubble wall; a set of 
timescales r oq , associated with different interactions that 
move the plasma toward chemical equilibrium or zero 
chiral charge; and tew, associated with the conversion of 
nieft into baryon number by the electroweak sphalcrons. 
As we show below, typically T cq <C r^is <C tew, and the 
dynamics of the first and second generation (s)fcrmions 
largely decouple from those of the third generation, which 
then become the dominant source of ne/s. 

We begin with the largest time scales. It is well known 
that the electroweak sphaleron time scale is tew ~ 
F^L - 10 5 /T, since T EW = §na 5 w T, with k ~ 20 
[H|, and aw the SU(2)i analog of the fine structure 
constant. The diffusion time depends on an effective 
diffusion constant for the plasma, D ~ 50/T [ill Il7l 
and the velocity of the advancing wall, v w ~ 0.05 ]Y2\ : 
r diff = D/vl~lO A /T^. 

The following reactions determine T cq : a) for third 
generation fermions, Higgs scalars, and their superpart- 
ners, Yukawa interactions associated with the decay, ab- 
sorption, and scattering of particles within the thermal 
plasma; b) strong sphalcrons that favor a relaxation of 
chiral charge to zero; c) supergauge processes involving 
spontaneous emission and absorption of gauginos, such as 
q + V «-» q. For example, the Yukawa-induced equilibra- 
tion time-scale for third generation left-handed quarks 
(q), right-handed stops (i), and Higgsinos (h), driven by 
the scattering q + t <-> h is numerically ~ 20/ YfT 
where Y t is the top Yukawa coupling for « 200 GeV 
and mi rts 100 GeV. Since Y t ~ 1, -C Tdiff, which 
in turn implies the approximation p q + fir — p^ = 
on Tdiff time scales. Similarly, the time scale for strong 
sphaleron-induced relaxation of the total chiral charge, 
N 5 = £j (2p qj - p Uj - p dj ), where qj, Uj, and dj denote 
the left-handed quark doublet and right-handed quark 
singlets of generation j, is T^ q ~ T" 1 ~ 300/T, since 
r ss = 6k' (8/3)a«T with a s being the strong coupling 
and k' ~ C(l) [HI- Hence, T| q <C Tdiff, leading to the 
condition N$ = on Tdiff time scales. 

Finally, when the masses of gauginos (V") are suffi- 
ciently light, supergauge processes can lead to chemical 
equilibration between SM particles and their superpart- 
ners, a situation we denote as "superequilibrium" defined 
mathematically as fi p = Hp = pp. This implies 

P=p + p={kp + kp)fi p -—=kpfi p —, (2) 

where k p are statistical weights that relate charge den- 
sity and chemical potential of particle species p: p = 



jp = kpfj.pT 2 /6. Even when gauginos become heavy, su- 
persymmetric Yukawa interactions can effectively bring 
about superequilibrium for third generation (s) fermions 
and Higgs(inos) in some regions of MSSM parameter 
space. The corresponding rates computed in [ltj Il7j in- 
dicate that typically T cq <C Tdiff for each supersymmet- 
ric three-body process. In the remainder of our analytic 
discussion, we will work in the superequilibrium regime 
where Eq. ([2]) holds. For example, we then have [ii = [it 
for (s)tops and [i~ h = [i h for Higgs(inos). 

Because the b quark and r lepton Yukawa couplings 
are small compared to Y t in the SM, they have been pre- 
viously neglected in EWB computations. In the MSSM, 
however, they can be significantly enhanced by tan /3. 
The lower bound on the mass of the lightest supersym- 
metric Higgs boson and electroweak precision data such 
as the anomalous magnetic moment of the muon favor 
tan/3 » 1. For tan/3 = 10, for example, Y^ T /Y 2 is 
100 times larger than in the SM, suggesting the possibil- 
ity of a relatively larger impact of the bottom and tau 
Yukawa interactions than in the SM [l3[ . After perform- 
ing a numerical scan over the relevant MSSM parameters, 
we find that the bottom (tau) Yukawa-induced equili- 
bration timescale is typically short compared to Tdiff for 
tan/3 > 5 (20) [TtJ- The Yukawa rates for the other gen- 
erations of fermions are typically too slow compared to 
the diffusion rate to have a significant impact. In what 
follows, we will work in the regime where t^*, t^?, and 
T^q are short compared to Tdiff and subsequently discuss 
possible departures from this domain. 

The timescale hierarchy has several physical conse- 
quences some of which have already been noted. First, 
by the time particles have diffused well ahead of the ad- 
vancing bubble, Yukawa-induced chemical equilibration 
and chiral charge relaxation have occurred. Second, for 
times < Tdiff, total baryon and lepton number are ap- 
proximately individually conserved, since their difference 
is always conserved and their sum is violated only the 
longer timescale tew- Focusing first on the Higgs scalars 
and third generation fermions, we observe that Yukawa- 
induced chemical equilibrium and supergauge equilib- 
rium implies to a good approximation that 

Llq+Vh-Vt = 0, flq-flh-Vb = 0, Hl-pLh-^T = (3) 

where I is the left handed lepton, the difference of the 
sign of [ih in Eqs. Q follows from conservation of hyper- 
charge. Adding the first two implies that 2[i q — p t — pi, = 
0, such that the third generation contribution to N$ van- 
ishes. As previously noted, on the time scale of Tdiff, we 
have A5 = 0. Moreover, because the first and second 
generation Yukawa couplings (diagonal or off-diagonal in 
gauge eigenstate basis) are tiny compared to the diagonal 
values for the third generation, there exist no significant 
interactions to generate non-vanishing first and second 
generation quark densities. Consequently, the first two 
generation densities approximately vanish for large tan (3. 

We now use the approximate conservation of baryon 
and lepton number on Tdiff time scale to relate the total 
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fermion plus sfcrmion densities to those for the Higgs plus 
Higgsinos. Since the first and second generation (s)quark 
densities vanish, the conservation of baryon number in 
terms of chemical potentials is fcQ/ig = —kpfJ-r — kpfJ-B- 
Using /ir + Ms = 2/iQ as implied by Eqs. (|3I2[) . we obtain 

k B - k T _ k Q k B - k T 
l l Q = ,. — I i. Mg -^Q = —- : —H 



C Q "I" K B ~r K T 



kn k Q + k B + k T 



(4) 

Similarly, applying lepton number conservation and 
Eq. ||2J) to third generation leptons implies 



L 



k H k L 



-H 



(5) 



where L (r) denotes the third generation left-handed 
(charged right-handed) lepton supermultiplet density. 
Using the approximate vanishing of first and second gen- 
eration fermion densities we obtain the total left-handed 
fermion density 



"loft 



kg k e 
m = j^-Q + —L 
kQ k L 



(0) 



kr 
k H 



H 



Eq. © is the key analytic result for the large tan /3 su- 
perequilibrium regime. It relates the non- vanishing Higgs 
supermultiplet density induced by CP-violating trans- 
port dynamics to m e ft via the statistical weights k p for 
third generation SM fcrmions, Higgs bosons, and their 
supcrpartners. The first and second terms on the right 
hand side correspond to the contributions from third left- 
handed generation quarks and leptons, respectively. 

The dependence on these contributions differs strik- 
ingly from what has appeared previously in the litera- 
ture: nioft ~ 5Q + 4T, with e.g., k B — kp k B — 9kp 
and kQ + k B + kp —> 9&Q + k B + 9kp in the numer- 
ator and denominator, respectively, of Eq. Q. These 
differences arise from (1) the contributions from the LH 
third generation (s)lepton density engendered when Y T is 
sufficiently large, and (2) the Y^-induccd bottom quark 
superequilibrium conditions of Eqs. ([3]). As discussed 
above, including the latter implies vanishing third gener- 
ation contribution to N§ and a corresponding suppression 
of the first and second generation quark densities. Using 
our full numerical solutions, we have verified that in the 
regime of small Y& and Y T , one recovers the previously 
identified dependence of ni e ft on Q and T . 

The presence of the statistical weights in Eq. ^ im- 
plies a strong dependence of Ub/s on the masses of the 
third generation squarks and sleptons, as one may ob- 
serve from the expression for the k p : 



k p (m p /T) _ c F ,i 



fcp(0) 



m p / 



dx- 

T (e x ± l) 2 



(7) 

where fc p (0) = 2g p (g p ) for Dirac (chiral) fcrmions and 
complex scalars with degeneracy g p , and cp( B ^ = 6(3) 



for fermions (bosons) and the + (— ) sign for fermions 
(bosons). For a fixed right handed (RH) stop mass, there 
exists a regime for « m t ~ in which k B rj Ht in Eq. j6]), 
corresponding to a nearly vanishing quark contribution 
to nieft- Indeed, this is the precise region in which the ad- 
vertised sign flip can occur depending on the magnitude 
of sparticle masses. Furthermore, if the quark contribu- 
tions are small, for sufficiently light tr, EWB may actu- 
ally be driven by the (s)lepton sector of the MSSM. This 
possibility does not arise under the previously studied as- 
sumptions of negligible Yf, tT and represents a qualitatively 
new class of supersymmetric baryogenesis scenario. 
Numerical solution and n B /s. Before turning to full nu- 
merical examples, we first isolate the essential physics 
with an approximate, analytic solution. In the regime for 
which Q and L are proportional to H as in Eqs. (|4l5p . it 
is convenient to combine the QBEs into a single equation 
for H. Eliminating all terms containing the fast Yukawa 
and strong sphaleron rates, we find 

-H Sf + Sf-Sf-Sf 
d x j^ = -T—+ H \ » , (8) 

k H 1 + K T + K L - K B 

where f = (T h + T T + T B + T r )/(1 + K T + K L - K B ) 
with Kp = H/P, the Tp being chiral relaxation trans- 
port coefficients that vanish outside the bubble Q, and 
where the denote the CP-violating source terms aris- 
ing from the scattering of supcrpartners p from the space- 
time dependent Higgs vevs. The diffusion ansatz allows 
one express the LHS of Eq. (|5J) in terms of the density H : 
dxjfj = H — DV 2 H, where D is the effective diffusion 
constant introduced earlier. 

Here, we will rely on the popular Higgsino CP vio- 
lating source (e.g. H, 0, Q) in order to illustrate the 
impact of the third generation sfermion masses and de- 
fer a study of possibly important S?^ and S?^ contri- 
butions. For simplicity, we assume a common relative 
phase M = Arg(/zM;6*) (i = 1,2), where b is the SUSY- 
brcaking Higgs mass parameter. As a benchmark, we 
work in the resonant regime and choose = \M2 \ = 200 
GeV and |M X | = 100 GeV, computing Sf and f from 
the expressions in Ref. 0] and a bubble wall profile given 
in [|| for tan (3 = 20. We consider a scenario with a light 
RH stop with mjj = 150 GeV as needed to obtain a 
strong, first order phase transition. All first and second 
generation sfermions and third generation LH sfcrmions 
have masses equal to one TeV. 1 

The results for ne/s as a function of the RH bottom 
squark mass, mg i , in units of the WMAP central value 
for sin^ = — 1 arc given in Fig.[TJ The solid and dashed 
curves give the results obtained with the complete nu- 
merical solution of the QBEs and the large tan/3 ana- 



The mass denotes the RH top squark mass at T = 0, after 
EWSB, while m t - is the finite T mass in the unbroken phase; 
similar notation applies for sbottom and stau as well. 
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FIG. 1: Baryon asymmetry in units of the WMAP central 
value as a function of RH bottom squark mass. The two solid 
(dashed) curves were obtained using the full numerical (large 
tan/3 analytic) solution. Light green (dark blue) curve corre- 
sponds to RH stau mass of 90 GeV (1 TeV). The dotted curve 
gives results neglecting bottom and tau Yukawa interactions. 

lytic solution embodied in Eqs. (|6|8|) . respectively, for 
two representative RH tau slepton mass: = 90 GeV 
(light green) and one TeV (dark blue). The dotted curve 
shows the numerical result obtained when bottom and 
tau Yukawa interactions are neglected. 

The impact of including realistic bottom and tau 
Yukawa interactions in the large tan/3, superequilibrium 
regime is striking, as seen in Fig. [T] Neglecting Yb yT gen- 
erally leads to a baryon asymmetry that is larger in mag- 
nitude, positive, and relatively insensitive to mg . The 
close agreement between the solid and dashed curves in- 
dicates our foregoing analysis captures the primary fea- 
tures of the dynamics in this regime. As expected, ub/s 
is largest in magnitude when is heavy and decreases 
as is decreased, reflecting the growing importance of 
sbottoms as they become light and the greater cancella- 
tion between kx and ks in Eq. (J5]). For heavy RH staus, 
the magnitude of the (s)lepton contribution to m e ft is 
small, and so ub/s vanishes when m~ h ~ mj and ks ~ fey. 



For light staus, the (s)lepton contribution becomes im- 
portant, and for either very heavy or very light sbottom, 
this contribution can change ti-q/s by a factor of two. 
(Since n^/n\ c it < 0, very light b and f lead to negative 
Ub/s, and for large m^ T , hb/s is positive.) 

Deviations from Eq. ^ can occur for small tan/3 or 
if superequilibrium is not valid. The region of tan/3 < 
10 corresponds to an interpolation between the upper 
dashed and the dotted curves in Fig. [IJ reflecting the 
larger third generation quark and smaller leptonic com- 
ponents to nieft- Although we have estimated that the 
off-diagonal CKM matrix elements are sufficiently small 
to justify neglect of flavor-changing 2— »2 scattering pro- 
cesses, such corrections should be kept in mind. We have 
also made (standard) assumptions regarding the trilinear 
scalar couplings which affect the decoupling of the light 
generations. Exploration of these and other effects will 
appear in forthcoming publications [l7| . 
Summary. During the upcoming era of more sensitive 
electric dipole moment searches and LHC studies, it is 
crucial to explore testable EWB scenarios such as MSSM 
baryogenesis [3| ■ We have shown in this context that the 
bottom Yukawa coupling is important even for moder- 
ate values of tan /3, which leads to dramatic changes in 
the basic physical EWB mechanism and the associated 
MSSM parameter space constraints. We also showed that 
the magnitudes of the sparticle masses can change the 
sign of the baryon asymmetry, and identified a new lep- 
ton driven supersymmetric baryogenesis scenario which 
does not involve RH (s)neutrinos. 
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